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We report the results of high-accuracy controlled laboratory measurements of the Stokes reflection matrix for 
suspensions of submicrometer-sized latex particles in water and compare them with the results of a numerically 
exact computer solution of the vector radiative transfer equation (VRTE). The quantitative performance of the VRTE 
is monitored by increasing the volume packing density of the latex particles from 2% to 10%. Our results indicate 
that the VRTE can be applied safely to random particulate media with packing densities up to ~2%. VRTE results 
for packing densities of the order of 5% should be taken with caution, whereas the polarized hidirectional reflec- 
tivity of suspensions with larger packing densities cannot be accurately predicted. We demonstrate that a simple 
modification of the phase matrix entering the VRTE based on the so-called static structure factor can be a promising 
remedy that deserves further examination. 
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Despite its extensive use in numerous areas of science 
and engineering for more than a century, the radiative 
transfer equation (RTE) had not had a solid physical 
foundation until quite recently. The derivation of the 
RTE for sparse discrete random media directly from the 
Maxwell equations [1,2] has finally made the RT theory 
(RTT) a legitimate branch of physical optics. Yet some of 
the assumptions used in this derivation are still formu- 
lated qualitatively rather than quantitatively. Perhaps 
the most important assumption is the limit of vanishing 
volume packing density of the particles populating the 
scattering medium: p -> 0. There has been some progress 
in quantifying this limit based on comparisons of p- 
independent RTT calculations with direct computer sol- 
utions of the Maxwell equations for random multiparticle 
groups with substantial packing densities [3—7] . How- 
ever, such comparisons are still impracticable for realis- 
tic media composed of extremely large numbers of 
particles, e.g., particle suspensions and powder surfaces. 
This makes controlled laboratory experiments an essen- 
tial tool for determining the maximal allowable packing 
density p max defined such that for p < p max the RTT pre- 
dictions can be expected to be quantitatively accurate. 

In this Letter, we report the results of controlled labo- 
ratory measurements of the Stokes reflection matrix for 
suspensions of microscopic latex particles in water and 
compare them with the results of numerically exact com- 
puter solutions of the vector RTE. By increasing the vol- 
ume packing density of the latex particles, we monitor 
the quantitative performance of the RTE and establish 
its range of applicability. We then discuss a simple 
“patch” that could be used to extend this range to larger 
packing densities. 

Note that the use of particle suspensions in this study 
was advantageous for two key reasons. First, it allowed 
us to vary the particle packing density in a controllable 
way (cf. [8,9]). Second, random Brownian motion of par- 
ticles during the measurement ensured the compliance of 


the scattering sample with the requirement of ergodicity 
necessary for the derivation of the RTE [1], 

The scattering geometry shown schematically in Fig. 1 
was designed to simplify the numerical solution of the 
RTE by making applicable the assumption of a laterally 
infinite plane-parallel scattering layer. Specifically, the 
scattering sample was an aqueous suspension of nearly 
monodisperse latex spheres filling a 25.4 mm diameter 
cylindrical cavity of 2 mm depth with absorptive side 
walls and bottom surface. The sample was illuminated 
by a parallel beam of light with a diameter exceeding that 
of the cavity, while the receiving optics sampled a small 
2 min diameter area near the sample’s center, as shown 
by the inset in Fig. 2. 

The nadir angle 8 0 of the incidence direction was kept 
constant at 20°. All scattering measurements were taken 
in the vertical plane through the incidence direction and 
were limited to the following zenith angles 6 of the re- 
flected light: -60°, -50°, -45°, -5°, 5°, 30°, 45°, and 60°; 
the sign convention is explained in Fig. 1. 

In the standard framework of plane-parallel scattering 
geometry, the polarized reflectivity of the sample is fully 
described by the 4x4 reflection matrix R((/. d 0 ) specify- 
ing the transformation of the Stokes parameters of the 
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incident beam into those of the reflected light [1]. The 
laboratory multidirectional spectro-polarimeter de- 
scribed in [ 10 ] was partially reconfigured to accommo- 
date an open, horizontally oriented hquid sample. This 
instrument allows one to measure with high accuracy 
the ratios R;j{0. 6 {) )/R n {0, 0 {] ) over a broad and quasi- 
continuous range of visible and near-infrared wave- 
lengths, where i,j = 1, ...,4. For the purposes of this 
study, the measurements were averaged over two 
10 nm wide intervals centered at the wavelengths 595 
and 695 nm. In agreement with theory, the elements of 
the upper-right and lower-left 2x2 submatrices of the re- 
flection matrix were found to be essentially zero. The 
residual deviations of these elements from the theoretical 


value zero were indicative of ~1% measurement uncer- 
tainties in the ratios Rij/R n . Therefore, the latter will 
be indicated by ± 1% error bars in the forthcoming plots 
comparing theory and experiment. 

According to the manufacturer (Bangs Laboratories, 
Inc.), the average radius of the latex microspheres was 
close to 350 nm, while their refractive index relative to 
that of water is known to be close to 1.196 at the two 
selected wavelengths. The measurements of the ratios 
Ry/R n were performed for three increasing values of 
the volume packing density: p = 2%, 5%, and 10%. 
Figure 2 shows a representative subset of the results thus 
obtained. Given the independent detailed knowledge 
of the scattering sample’s macro and microphysical 
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Fig. 2. Measurements (filled circles), RTT computations (dashed curves), and RTT+SSF computations (solid curves) of the ratios 
of the elements of the Stokes reflection matrix R for aqueous suspensions of latex particles with packing densities p = 2% 
(red color), 5% (yellow color), and 10% (blue color). The two gray curves show the best fit of RTT + SSF theoretical results 
for p = 15% to the laboratory data for p = 10%. The inset is a photograph of the scattering sample in which the red dot shows 
the size and approximate location of the area sampled by the receiving optics. 
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parameters, this dataset constituted an excellent quanti- 
tative test of the RTT. 

To solve the vector RTE for the plane-parallel geom- 
etry, we used the well-known adding/doubling method 
[1,11] in which the reflection and transmission matrices 
of the water-air interface were incorporated according 
to the Kirchhoff approximation [12], otherwise known 
as the ray-optics approximation [13]. Based on extensive 
calculations, we found that below a certain nontrivial 
surface roughness threshold, the cumulative reflection 
matrix R(YL 0 0 ) of the entire system becomes essentially 
independent of the surface roughness scale except at ze- 
nith angles 9 extremely close to the specular reflection 
direction. This implies that microscopic undulations of 
the water surface that could be caused by air pressure 
fluctuations and/or mechanical vibrations of the experi- 
mental setup had negligible effect on our measurements 
and calculations for the specific sampled scattering 
directions. The single-scattering quantities entering the 
RTE were computed according to the standard Lorenz- 
Mie theory [14], 

Instead of using the manufacturer-specified character- 
istics of the latex-particle suspensions, we considered 
them unknown and “retrieved” them by looking for the 
best fit of theoretical RTE results to the laboratory data 
at both selected wavelengths. We found that the most dif- 
ficult task was to fit the measured ratio R 44 /R n . This ne- 
cessitated, for example, the use of an extremely narrow 
size distribution of latex particles. Nevertheless, as Fig. 2 
demonstrates, a near-perfect fit was obtained in the case 
of the 2% suspension for the following values of the 
model parameters: effective radius of the size distribu- 
tion 350 nm, effective variance of the size distribution 

0.001, refractive index 1.196, and optical thickness 128 
at 595 nm. All these values are in excellent quantitative 
agreement with the independently available information. 
Importantly, no satisfactory fit was found for effective 
radii and refractive indices deviating by more than 
±15 nm and ±0.04 from the corresponding best-fit 
values. 

Obtaining a good theoretical fit to the measurement 
results for the 5% suspension was also possible and re- 
sulted in the same values of the microphysical latex- 
particle parameters (see the two R 44 /R u panels in Fig. 2). 
However, the corresponding optical thickness value 485 
at 595 nm was significantly greater than 2.5 times that for 
the 2% suspension. Fitting the laboratory data for the 10% 
suspension turned out to be impossible even when the 
model optical thickness was increased to infinity. 

To improve the fits for the 5% and 10% suspensions, we 
modified the phase matrix entering the classical vector 
RTE by applying a simple phenomenological patch based 
on the so-called static structure factor (SSF) [15,16], This 
phenomenological recipe is intended to partially account 
for mutual correlations of particle positions in a densely 
packed discrete random medium and serves to make p an 
explicit model parameter. The SSF was computed using 
the Percus-Yevick approximation for impenetrable 
monodisperse spheres [17]. 

Figure 2 shows that despite its technical simplicity, the 
introduction of the SSF appears to be a promising mod- 
eling tool. Although the fit remained virtually the same 
for the 2% and 5% suspensions, the scaling of the optical 


thickness became more realistic (115 versus 280). The 
solid blue R 44 /R n curve at 595 nm is now completely 
within the measurement error bars for an optical thick- 
ness of 550, and the corresponding fit at 695 nm was sub- 
stantially improved. 

Interestingly, a nearly perfect fit to the laboratory data 
for the suspension with the 10% packing density can be 
obtained with RTT ± SSF computations for p = 15%, as 
the two gray curves in Fig. 2 demonstrate. This result sug- 
gests that the RTT ± SSF patch may necessitate values of 
the formal packing-density parameter somewhat larger 
than the actual ones to yield an improved fit to the results 
of bidirectional photo-polarimetric measurements. 

In summary, our results further corroborate the meth- 
odology used in the microphysical derivation of the vec- 
tor RTE from the macroscopic Maxwell equations [1,2] 
and suggest that the RTE can be applied safely to random 
particulate media with packing densities up to ~2%. RTE 
results for packing densities of the order of 5% should be 
taken with caution, whereas larger packing densities can- 
not be accurately handled and definitely require the use 
of alternative modeling tools [18]. In this regard, the sim- 
ple modification of the phase matrix based on the SSF 
shows certain promise and should be further analyzed. 
Given the widespread (and often uncritical) use of the 
RTE in many science and engineering disciplines (e.g., 
[19-25] and references therein), the results of our study 
should provide important guidance in the process of 
identifying appropriate modeling tools and deciding 
whether the RTE can even be used depending on the spe- 
cific practical problem in hand. 
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